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ABSTRACT 



Drag of free-falling spheres in water and in dilute aqueous 

solutions of poly (ethylene oxide) was measured by ejecting the 

spheres near terminal velocity in a fluid-filled tank. Polymer 

additive concentrations tested were 1000 wppm, 200 wppm, and 100 

wppm. The results agree with previous investigations, available 

only for subcritical Reynolds nunibers, in that the drag coefficient 

decreases with increasing water Reynolds number and decreasing 

polymer concentration. It is shown that in the polymer solutions 

abrupt decrease in drag coefficient occurs at about the same 

5 

Reynolds number (3 x 10 ) as the critical value in water. For 
Reynolds numbers greater than the critical value, the drag was 
still less than in water and increased slightly with increasing 

5 

polymer concentration. At Reynolds nuttibers between 7 and 8 x 10 , 
the drag coefficients for all polymer concentrations, including 
water, abnptly increased to nearly a canmon value at the upper limit 
of this investigation. 
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1 . Introduction 



1.1 General . Drag reduction (or friction reduction) by use of 
polyiier additives in a flowing liquid has only recently begun to 
be developed. In 1944, Vanoni (61) observed that in a flune 
suspended sediment tended to ^preciably reduce the turbulent 
transfer of momentum and hence the resistance to flow. This effect, 
which seerred to depend vpon sediment concentration, allowed the 
sediment laden water to flow more rapidly than a coitparable clear 
water flow. Then, in 1948, Toms (51) reported correlation between 
concentration and increased flow-rates (at constant pressure gradient) 
for turbulent flow of a high molecular-weight polymer solution in 
a circular pipe. Since this was apparently the earliest positive 
correlation between friction reduction and the concentration of 
polyiter additive, the effect was termed the "Toms' phenomenon" by 
various authors. After Tom's report, although considerable effort 
was given to the study of the rheological properties of non-Newtonian 
fluids, little interest was evidenced in literature concerning the 
friction reduction properties of these substances until the late 1950 ' s 
when the possibilities of practical uses for the Tom's phenomena began 
to be realized. Petroleum coitpanies becaite interested when it was 
demonstrated that by introducing minute quantities of certain high 
molecular-wei^t polymers into pipelines the power required and thus 
the cost to punp petroleum or other fluids could be considerably 
reduced (7, 58). Other chemical process industries became interested 
in friction-reducing additives for the same reason. 
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1.2 Waterborne Vehicles. Another obvious potential for friction 



reduction was in applications to waterborne vehicles and the works 
of many researchers are reported in shipbuilding journals and 
reports of agencies connected with naval engineering. In order 
to study friction reduction on water vehicles, the problem was 
sometimes regarded as friction of a flat plate plus drag due to the 
wake generated as the vehicle moved through the water. Granville (17) , 
assuming that practical drag reducing substances were non-Newtonian, 
applied ideas that had been developed by Metzner and his associates 
(32) for non-N©A?tonian additives in pipe flow to flat plates. He 
predicted a decided drag reduction for flat plates in a flowing non- 
Newtonian power-law fluid. Vogel and Patterson (62) reported drag 
reduction resulting from injecting poly (ethylene oxide) , ( POLYO)0 , 
solutions into the boundary layer of a three-dimensional streamlined 
model. They correlated drag reduction with polymer concentration 
and observed that the drag reduction increased with increasing 
molecular weight of the linear polymer used. Kowalski (25) observed 
significant drag reduction on ship models and flat plates towed in 
a tank containing very dilute solutions of POLYOX. Although he 
referred to this fluid as being non-Newtonian conventional viscometry 
indicates that even at high rates of shear these dilute solutions of 
polyox are Newtonian. A further demonstration of drag reduction by 
"non-Newtonian" additives for water vehicles was given by Thurston 
and Jones (50,18) who tested a torpedo shaped vehicle having a 
soluble coating of high molecular weight polymer on its nose section. 
Merrill, Smith, and Chung (31) reported drag reduction on both a flat 
plate and a torpedo shaped body for low concentrations of polyox in 
water. However, drag on -the "torpedo" increased for high polymer 
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concentrations. This was explained by assuming that the boundary 
layer had been laminarized so that early separation occurred and 
thus a larger diameter wake led to increased pressure drag. 

Although, in principle, drag reduction for water-borne craft 
can be accoitplished, utilization of the principle on a full scale 
may not be practical unless the efficierjcy of the process is such as 
to allow practical quantities of polymer to be eitployed. Kowalski 
(25) estimated that to reduce drag on a typical 450 foot merchant ship 
at a speed of 18 knots would require injecting on the order of 
13000 Ib/hr. He pointed out the obvious fact that continuous application 
of this additive would not likely be at all practical. However, the 
use of additives for short bursts of speed in large vessels or continuously 
in much smaller vehicles might be more feasible. These exarrples of work 
in the field of drag reducing additives indicate that airple motivation 
is present for continuing study of these si±>stances. By gaining more 
insight into the nature of turbulent friction reduction not only will 
more practical utilization of the principle be realized but it is likely 
that greater understanding of turbulence itself may result. 

1.3 Fluid Dynamic Drag and Drag Reduction . The drag force or resistance 
to motion of a totally submerged body moving through a fluid is usually 
thought of as consisting of resistance due to skin friction plus profile 
drag due to a negative pressure developed in the wake behind the body. 

Skin friction and profile drag, taken together, are often referred to 
as viscous drag. In the case of bodies moving in water and which are 
not fully submerged an additional drag or resistance component known as 
wave drag affects the motion of the boc^. Wave drag results from surface 
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wave effects and surface spray effects. Since men began using ships 
and various water vehicles, attention has been given bo methods of 
reducing the resistance to notion of the vehicles through water. Until 
recently the efforts have been largely confined to wave and profile drag 
by reducing frontal area or to reducing skin friction by smoothing 
surfaces ej<posed to the fluid and by controlling the fluid flow around 
the body so that the boundary layer remained laminar. As mentioned above, 
the possibility of using polymer additives to reduce the viscous drag 
has generated considerable interest in the investigation of fluids 
containing these additives. 

The analysis of fluid flow past various simple geometrical sh^jes 
has often provided insight into solution of problems of flow past more 
conplicated but more practical shapes. The resistance to flow past a 
flat plate or resistance to flow in pipes is almost entirely due bo 
skin friction. There is no profile drag due to a wake nor is there a 
boundary layer separation. On the other hand, except for flows with 
very low Reynolds number, profile drag is predominant for flews past 
spheres and past cylinders oriented normal to the direction of flow. 

Thus, by a suitable choice of mcdel shapes various contributions to 
the resistance to fluid flow can, to some degree, be "isolated". 

Although this investigation is primarily concerned with drag 
reduction for spheres, a brief summary of those fluids flow char- 
acteristics affecting drag for flows in various geometries is presented 
in order that drag or friction reduction may be viewed in broad perspective. 
In particular, much emphasis is given to pipe flow, largely because of 
the abundance of work reported for friction reduction in pipe flow. 

This, in turn, is due to the relative ease in performing pipe flow 
measurements and the difficulty 
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of cxjnducting ej<periraents with other shapes. In addition, topics 
discussed in the sunmary will enhance understanding and interpretation 
of sphere drag reduction. 

1.3.1. Rheology . The relationship ^called the constitutive 
equation , between flow paraitieters and the shear stress exerted by the 
fluid upon a bounding surface or between adjacent fluid elements is 
an essential element in an analysis of viscous drag. For those fluids 
which are called Newtonian and can be considered inooirpressible this 
relationship for a one-dimensional flow is: 

X = y(du/dy) 

viiere x is the shear stress, u the constant shear viscosity of the 
fluid, u the oonponent of velocity normal to y, the direction of 
displacement between the fluid element and the bounding surface or 
bebveen the adjacent fluid elements. A plot of shear stress with 
rate of shear, t versus du/dy, for Newtonian fluids is linear. Those 
fluids which do not have this linear flew curve are collectively 
called "non-Newtonian". Flow curves for some fluids which are 
encountered among various time-independent fluids are shown in 
Fig. 1.1. See Skelland (48). A Bingham plastic exhibits a constant 
fluid viscosil^ but only after a "break-out" shear is applied. The 
pseudoplastic flow curve characterizes most non-Newtonian fluids. 

These fluids are often called shear thinning, which seems appropriate 
ipon noting the shape of the flow curve. A few fluids, semetimes 
termed "shear thickening", are characterized by the "dilatant" fluid 
flow curve type. Many fluids, over a restricted range of shear rates 
can often be characterized by the so-called power-law or Ostwald- 
de Waele model: 

X = K (du/dy)’^ 



where n and K are ej^rojtentally determined oanstants. Obviously 
this model is valid only for those values of shear rates vhich 
produce linear log t versus log du/dy plots. Newtonian fluids 
represent the special case of a power- law fluid with n = 1 and 
K the shear viscosity. Shear thinning fluids have 0 < n < 1 and 
K can be thought of as an apparent viscosity although, since this 
is eitpirical, no physical significance can be inferred. Another 
inportant groip of non-Newtonian fluids exhibit time dependence in 
that they have both viscous properties and elastic properties. 

When subjected to stress these fluids flow, but when the stress is 
removed they tend to partially recover. For these fluids the 
constitutive equation given by Skelland (48) and attributed to 
Frohlich and Sack (11) and Oldroyd et. al. (35) is 

t + A^(d./dt) = y t (dy/dt) + A2(d^Y/dt^)] 

where and are experimentally determined relaxation times, y 

* 

IS the strain, and u is an experimentally determined viscosity 
which at very low rates of shear is the Newtonian viscosity. It is 
noted that for a purely viscous fluid, i.e., a flmd having no 
significant elastic properties, this equation reduces to that given 
above for the Newtonian case since a^ = A 2 = 0 and for low strains 
dy/dt = du/<^. For high rates of shear, a constitutive equation 
containing higher order terms and five additional experimentally 
determined constants has ixeen derived by Oldroyd et.al. (35) . 

The drag reduction phenomenon produced in fluids containing 
even trace quantities of polymer additives are not readily explained 
by accepted descriptions of Newtonian flew. Thus, these fluids are 
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classified as non-NewtxDnian, although there is evidence (9) that, 
in terms of conventional viscometry, drag reduction occurs in sorte 
fluids which appear Newtonian. 

1.3.2. Pipe flow . For flow in circular pipes a friction 
coefficient is defined as: 

A = (-dp/dx) D/(jpV^) 

where dp/dx is the drop in pressure per unit length of pipe, D the 
pipe diameter, p the fluid density, and V the mean fluid velocity 
in the pipe. If a Reynolds number is defined as Re = VD/v where v 
is the Newtonian kinematic viscosity of the fluid then it has been 
shown, both by solution of the Navier-Stokes equations and by 
ej<periment, that for laminar flow of a Newtonian fluid: 

A = 64/Re 

This IS often called Hagen-Poiseuille flow. This ej^ression for pipe 
friction coefficient generally has been obseirved to be valid for 
Reynolds numbers less than about 2100 although the exact transistion 
Reynolds number depends ipon the degree of success achieved in injecting 
the fluid into test sections without generating flow disturbances 
leading to turbulence. 

For fully developed turbulent flew in pipes, no general and exact 
solution of the Navier-Stokes equations is available; however, semi- 
eitpirical expressions have been developed. Prandtl's mixing length 
hypothesis and his assumption of a simple form for the mixing length, 

1 = !Cy led to a universal velocity distribution law for the flow of 
Newtonian fluid characterized by large Reynolds members: 

(u/v*) = (l/K)ln(yVj^/v) - B 
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where u is the tune average fluid velocity at a point, v^ the 

1/2 

friction velocity = ('’^/p) / the tangential wall stress, 

y the distance from the pipe wall locating the point at which the 
velocity is u, and < and B constants of integration which must be 
determined experimentally. Prandtl's universal law of friction 
for srtooth pipes can be obtained from the universal velocity 
distribution law given above. 

= 2.0 log (Re - 0.8 

The constants in this ej^ression result from setting < = 0.4 and 
B = 5.5 in the universal velocity distribution law. These values 
produce e>p>ressions vhich best agree with velocity profiles and 
friction coefficients observed for the flow of water in smooth 
pipes. 

Dodge and Metzner (8) developed enpirical friction-coefficient 
ejpiressions for turbulent pipe flow for fluids vhose constituitive 
equation could be approximated by the power- law model. This enpirical 
e>pression was: 

l/fl/2 ^ log (Re'f^"”'^^) - 0.40/(n’)^*^ 

where n' is the flow behavior index, f is the Fanning friction 

n * * n * “J. 

factor if = a/4) , and Re' = (D v c)/(g k'8 ) is a generalized 

Reynolds number. The flow behavior index is determined experimentally 
from: 

n' = d(log(DAP/4L))/d(log(8V/D)) . 

The constants are further related by: 

= DAP/4L = K' (8V/D)’^' 
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where AP/L is pressure drop in a length, L, in the pipe. For pipe 
flow these constants can be related to the power- law constants by: 



Metzner and his associates (32) had previously shown that for laminar 

pipe flow, f = 16/Pe' or X = 64/Re'. Although the pcwer-law model has 

often been useful, many fluids of interest as drag reducers are not 

well described by this general egression (8,45). In addition, those 

fluids having significant elastic properties would not be expected to 

behave according to this description. 

As has been mentioned previously, many fluids in which drag 

reduction occurs are essentially Newtonian. Some researchers have 

attempted to determine new constants for the Prandtl friction law 

for various concentrations of additives in fluids (16) or have 

1/2 

correlated friction coefficients using, as parameters, 1/f and 



friction law of POLYOX WSR 301, Goren and Norbury (16) determined 
mixing length values for various additive concentrations and distances 
from the pipe wall. Many workers have suggested the likelihood (16,21, 
39) that the drag reduction observed for turbulent pipe flow results 
from viscoelastic effects, although the elastic effects are difficult 
or inpossible to observe by the usual methods. 

One of the problems associated with these investigations was that 
some polymer additives degraded in turbulent flow thereby altering the 
drag reduction properties. Seme of the fluids were clearly viscoelastic 
when containing polymer additives in excess of about 0.1% (1,000 weight 
parts per million, vppm) but the concentrations enployed in many of the 
investigations were on the order of only a few wppm. Because of these 



K' = ((3n + l)/4n)^ 




Incidental to determining empirical constants for the 
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problems and because of the apparent Newtonian character of many 
of the drag reducing fluids, the explanation of drag reduction by 
viscoelasticity has not been universally accepted. 

In an analysis of the effect of polymer additives on turbulent 
wall friction F. M. White (66) attributed to Meyer the observation 
that the constant, B, in the universal velocity distribution law 
for pipe flew, changes upon addition of polymer if the friction 

ic 

velocity exceeds sane threshold value, v^ . Meyer gave an ^proxi- 

* ★ 

mation for B as: B ^ 5.5 + aln(v /v ) 

where a is a dimensionless constant dependent tpon the type of additive 
and the concentration. The constant < is apparently unaffected by 

polymer additives. White gives a for POLYOX WSR 301 solutions: 

^ 0.57 
a = 3c 

where c is polymer concentration in vppm. Some of the typical friction 
coefficients observed for various fluid types flowing in pipes assumed 
to be smooth are given in Figs. 1.2, 1.3, and 1.4. The solid lines 
denote curves for pure solvent and the dotted lines represent behavior 
of the fluid with additives. In each case, it is significant that in 
the laminar flow region the additive does not change the fluid friction. 
Hershey and Zakin (21) noted that a fluid containing drag reducing 
additives remained in laminar flow until a greater Reynolds number had 
been achieved than for the pure solvent. For fully turbulent pipe flow 
a second kind of drag reduction takes place. In figures 1.3 and 1.4, 

A, B, and C represent different concentrations. The primed letters of 
figure 1.3 denote flow in smaller diameter pipe than that for the unprimed 
letters. The "critical" Reynolds number for drag reduction in turbulent 
flow has been shown to be dependent upon pipe diameter (9,21,39) and 
independent of polymer concentration (16) . Ram attributed this to a 
relaxation within the polymer chain. 
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1.3.3. Flow past flat plates . In the case of the t>?o-dirinensional 
boundaiy layer flow along a flat plate at zero incic3ence the friction 
coefficient is defined as: 

C = /i/2 oU ^ 



where is the fluid velocity far fron the plate. The velocity profile 
for smooth flat plates can be written in a logarithmic form similar 
to that given for turbulent pipe flow but it is usually more convenient 
to express flat plate resistance coefficients enpirically such as 
the Prandtl-Schlicting skin friction formula: 

= 0.455/ (log - A/R^ 

where is Reynolds number based on the length, L, of the plate; 

=! U^L/v. For laminar boundary layer flow along a flat plate, the 
Blasius formula, =1328/Rj^^^ ^ applies. Granville (17) developed 
friction coefficients for flat plates in power-law fluids. For 
laminar flow: 

= 2(0.824)^'^^’^’^^^ [(n+l)^'*'^^^/( 2 n+l) (3n+2) 

vdiere n is numerically equal to the value determined in pipe flew ej^riments 
Rj^' is the generalized Reynolds number based on length of the plate. 



For turbulent flow: 






Cg = pV'log R]-_ 



where P and Q are eirpirically obtained. Van Wijngaarden (60) as well 
as Acrivos et. al, il) following Granville, presented an analysis of 
laminar ix)undaty layers tor pewer-iaw fluids. F. M. t'hite (66) used 



23 



data from pipe flow to predict friction coefficients for turbulent 
flow along flat plates. He observed that, as shown in Fig. 1.5, 
reduction in local skin-fraction coefficient would be ej^ected in 
laminar flow as well as in turbulent flow, however, contrary to 
observation in pipe flow, the turbulent friction reduction would 
decrease with increasing solvent Reynolds nunber and with decreasing 
a, a paraireter determined frcm pipe flow data. At some point no drag 
reduction would occur for fiorther Reynolds number increase. 

1.3.4. Torque on Rotating Disks . A circular disc rotating 
about Its axis normal to the surface is a convenient device for 
investigating drag reduction. The torque on the disc is entirely 
due bo skin friction and this torque can be easily measured. For 
a rotating disc in a stationary fluid, a moment coefficient is defined 
as : 

2 5 

= (torque on both wetted sides of disk)/ (l/2pu; r ) 

with 0 ) the angular velocity and r the disk radius. Schlicting (43) 
gives for Newtonian laminar flow: 

C = 3.87/Re^'^^ 
m 

2 

where the Reynolds nunber. Re = r cu/v. For turbulent flow Goldstein 
(15) developed a moment coefficient for rotating disks of the same 
form as the universal resistance law for pipes: 

1/C ^ ^ ^ jpg Q 1/2 ^ 

m ^ m 

Hoyt and Fabula (23) found that, within the scope of their investigation, 
increasing the concentration of poly (ethylene oxide) up to 621 wppm 
produced increasing drag reduction. The trend, as shown in Fig. 1.6 
seems to be that turbulence is suppressed, moving the flow conditions 
toward laminar Newtonian .flow. 
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1.3.5. Flow Past Bluff Bcdies . Drag exerted by a fluid 



flowing past a bluff body is more difficult to analyze than 
friction in pipes or on flat plates. Therefore, Analytical 
ej^ressions for cylinder and sphere drag coefficients such as 
those by Lamb, Stokes, and Oseen have been developed only for very 
low Reynolds nunibers. Numerical solutions such as that for 
spheres given by Jensen (24) are available for Reynolds numbers 
vp to about 40. Beyond this, general analysis is not very practical. 

In pipe flow and flow past flat plates only skin friction was involved 
in the resistance to fluid motion. On the other, flows past bluff 
bodies at high Reynolds numbers are characterized by vortex formation 
and shedding accompanied by boundary layer separation. This leads 
to formation of a wake and a negative pressure v^ich acts as a 
retarding force, termed form or pressure drag. In addition, as 
Reynolds number is increased the boundary layer changes from laminar 
to turbulent. Both cylinders and spheres exhibit sudden decreases 
in drag coefficient at a Reynolds number referred to as the critical 
Reynolds nuirber. Drag coefficients are defined as: 

2 

= (drag force/unit length) /(l/2pU d) for cylinders, 

2 2 

Cj^ = (drag forceJ/(l/2pU ) (nd /4) for spheres, 

and Reynolds numbers are Re = Ud/ (for both cylinders and spheres) , 
vbere d is diameter and U is free stream fluid velocity. In both 
cases the critical Reynolds number occurs around 3 x 10^ the exact value 
depending upon surface roughness and turbulence in the fluid. The 
curves shewn in Figs. 1.7 and 1.8 are drawn to indicate most commonly 
encountered critical Reynolds numbers and largely represent the work 
of Weiselsberger as presented in current literature such as that of 
Schlicting (43) and Hoemer (22) . 
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For very low Reynolc3s numbers, skin friction is a significant 
carponent of drag for bluff bodies. For spheres , Stokes ' solution 
shows that two-thirds of the drag at low Reynolds nun±>ers is due 

to skin friction. Goldstein (15) gives an ejqsression for cylinder 

1/2 4 

skin friction as = 4/Re which is valid for 30 ^ Re < 10 and 
DF 

at Re = 10 Goldstein says Thom determined that skin friction accounts 
for about 43% of the drag. As Reynolds numbers are increased and 
wake size increases, pressure drag becomes the significant contributor 
to drag of bluff bodies. Upon further increase of Reynolds number, 
the boundary layer becomes turbulent, the point of separation of 
the laminar boundary layer from the bocfy surface is shifted abrtptly 
toward the rear, and drag is reduced by the resulting reduction in 
wake size. In the case of cylinders, Roshko suggested that this 
shift rearward is actually a reattachment of the boundary layer to 
the surface forming a bubble vdiich he called a separation bubble 
followed by a second separation further toward the rear. Sanders 
(42) suggests the possibility of the same mschanism working for 
spheres also. 

Because of the ability to assume a two-dimensional boundary 
layer /Cylinders, at first, appear simpler for investigation than 
spheres, but complications arise because of end effects and the 
periodic shedding of vortices which give rise to oscillatory forces 
transverse to the direction of flowing fluid. Consequently, an 
infinite length-to-diameter ratio (L/D) must be assumed to obtain a 
generalized plot of cylinder drag coefficients with Reynolds number. 

If one works with cylinders whose L/D is small he must either apply 
some correction factor or state the value of L/D for which his results 
are valid. 
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The degree of difficulty in handling Nev/tonian flows about 
bluff bodies suggests that to further conplicate the problem by 
disallowing the Newtonian assunption for a constituitive equation 
would bring about insurmountable analysis difficulties and, to date, 
this has been the case. Investigations of non-Newtonian flow 
about cylinders are represented by a few applications of two-dirrensional 
boundary layer solutions such as that by Acrivos, et al (1) for power- 
law fluids and a few experimental investigations of fluids containing 
drag reducing additives such as that by McClanahan and Ridgely (27) . 

The latter authors reported drag enhancement for Re less than about 

4 

10 and drag reduction for greater Reynolds numbers. Although the 
scope of their work did not extend to the critical Re, the trend of 
their results as sketched in Fig. 1.9 suggested that, in fluids 
containing drag reducing additives, there might be no "drag crisis", 
as the sudden drop in drag coefficient at the critical Reynolds nurrber 
is sometimes called. 

Spheres in non-Newtonian fluids have enjoyed slightly more attention 
than have cylinders. For very low Reynolds numbers, Skelland (48) gives 
for power- law fluids: 

= 24 X/(d^^“^p/K) 

where X is a function of n and is determined experimentally. This 
expression is not very satisfying in that many variations in X have 
been found by various investigators. Slattery and Bird (49) presented 
another analysis for sphere drag coefficients at very low Reynolds 
number for fluids which are characterized by the Ellis Model: 

I = (du/dy)/(A + 



21 



where A, B, and a are constants to be evaluated ertpirically . 

Bizzell and Slattery ^4) calculated boundary layer separation 
points on spheres as a function of n, the power-law model ej<ponent. 
Turian (57) applied the Ellis model to one to two percent aqueous 
solutions of poly (ethylene oxide), and hydroi^ethyl cellulose, in 
an investigation of spheres in creeping flow. These theoretical 
results provide some insight into behavior of non-Newtonian fluids 
but they do not allow prediction of the behavior of flmds containing 
drag reducing additives for flows characterized by fully developed 
wakes and separated boundary layers. 

Drag reduction has been reported by various researchers as well 
as drag enhancement in a few cases for spheres in fluids with polymer 
additives. These results have been summarized by Sanders (42) and 
the general trend is shown in Fig. 1.10. A. White (64) performed a 
particularly interesting experiment by roughening the surface of a 
sphere so as to artifically "trip" the boundary layer into turbulence 
and thereby produce minimum pressure drag in water . Dropping the 
sphere in a very dilute aqueous solution of poly (ethylene oxide) he 
observed a drag incnrease. His eiqplanation was that the polyirer 
additive had suppressed the IxDundary layer turbulence and caused the 
pressure drag to return toward the value of drag obtained prior to 
roughening the sphere surface. 

1.4 Drag Reducing Polymers 

1.4.1. General . In the foregoing summary of fluid drag it 
was apparent that most investigators had failed to correlate drag 
reduction with non-Newtonian behavior. In seme cases polymer solutes 
were sufficiently concentrated that the power “law model or Ellis model 
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cx)uld be ^plied with slight success but, in general, drag 
reduction had to be attributed to sorething other than just a 
variation in rheological properties. It has been suggested by Hoyt 
and Fabula (23) , Hershey (21) , and others, that a possible explanation 
might be that viscoelastic effects in even apparently Newtonian fluids 
are responsible for turbulence suppression and consequent reduction 
of skin friction of pipe flow, flow past flat platen and rotating disks. 

In the case of bluff bodies the viscoelastic effects would be assurted 
to reduce turbuleiice in the wake and del.^ separation of the laminar 
boundary layer causing reduced pressure drag. 

A wide variety of polyirers have been used as drag reducers. Among 
the most coirnionly used in water are poly leriiylene oxide) (PEO) , guar 
gum, polyacrylamide (SEPAEAN) , carbojqmnethyl cellulose (CMC) , and 
hydroxyethyi cellulose (HEC) . Polymers used in other fluids are 
poly (methylmethacrylate' (PMMA) in toluene and in monochlorobenzene 
and poly isobutylene (PIB) in kerosene and in cyclohexane. Hoyt and 
Fabula (23) reported results of drag reduction tests on twenty-eight 
different polymer additives and five natural si±)stances. Ihose 
producing the most dramatic drag reduction were PEO (WSR-205) , guar 
gum, and SEPARAN. For 70% torque reductions on a rotating disc at 

5 

Re = 1.3 X 10 , concentrations in wppm required were 250 for WSR-205, 

500 for guar gum, and 100 for SEPARAN AP-30. In other tests of the 
same investigation they observed that other forms of PEO having higher 
molecular weights, WSR-301 and COAGULANT, gave even more striking results. 
In concentrations in excess of about 0.1% many drag reducing additive 
polymers cause the resulting fluids to exhibit measurable viscoelasticity, 
and a significant portion of the wDrk done with drag reducers has been 
at these concentrations. Of special note, however, is that PEO and guar 
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gum, as well as several of the others, yield greater drag reduction 
when present in concentrations on the order of 100 wppm, and less, 
than when concentrations are such that viscoelasticity can be 
measured. For all geometries it has usually been reported that 
drag reduction increased with concentration to some maximum drag 
reduction, and then decreased with further increase in concentration. 

It has been generally observed that drag reducers have very 
long-chain molecular structures and have very large rtolecular 
weights, on the order of 10^. Quite appropriately they are often 
referred to as itecromolecules . High moleculcu: -weight polymers having 
chain structures with many appendages or highly coiled chains, generally 
have not produced as dramatic drag reductions as have those whose 
structures are simple (38) . The soli±ulity of the additive also 
affected the drag reduction (23) . It must be recognized that tem- 
perature and presence of inorganic salts often dramatically affect 
chain coiling v^ich in turn affects polymer solubility. Another, 
sometimes annoying, prcpetty of sore high molecular-weight polymers 
is their susceptibility to degradation. It has been found that many 
polymers ^e reduced in size by chain scission brought about by 
oxidation, turbulence, droplet formation, and sound waves. This 
degradation, if it occurs in the course of a drag reduction experiment/ 
often makes interpretation of the results difficult, especially if one 
IS attempting to correlate molecular weight with drag reduction. 
Determination of the degree of degradation during an experiment is not 
as easy task. One must first solve the problem of taking a representative 
sample and then determine molecular weight by seme means such as by 
determining intrinsic viscosity. The degraded polymer's molecular weight 
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IS then calculated from the widely accepted enpirical relationship: 



[r.] = kM 



,0.78 



where [r] is the intrinsic viscosity, k is a teirperature dependent 
parameter, and M is the molecular weight. The intrinsic viscosity 
IS defined as: 



Ixm 

C'O 



(('- - 






where n = u, the shear viscosity of the solution, is the solvent 

o 

viscosity, and c is the solute concentration expressed in grams per 



100 ml. 

1.4.2. Poly (ethylene oxide) . Poly » ethylene oxide), also 
called polyoxyethylene and polyoxiraie, is an addition polymer 
of ethylene oxide and has molecular weights ranging between 100,000 
and 10 million. Lower molecular weight ethylene oxide polymers 
are usually kncMn as polyethylene glycols. Union Carbide Chemicals 
Conpany presently is the sole manufacturer of PBD and the polymer 
IS marketed as "POLYOX" in several molecular weight ranges, some 
of v\hich are listed in Table I. Powell and Bailey (36j note that 
PED has several properties vihich vary from those expected of glycols 
and these apparently result from the extremely long molecular structure. 
The resin is water soli±)le at room tenperature but has an inverse 
tenperature solubility relationship; i.e., as tenperature increases, 
solubility decreases. Presence of various salts such as sodium 
phosphate and sodium carbonate can significantly lower the solubility 
of PBO in water. It is soluble in aromatic solvents but insoluble 
in aliphatic hydrocarbons and glycerin. Aqueous solutions of PBO 
at concentrations of 1% and greater are pseudoplastic and definitely 
are not thixotropic. Solutions with concentrations on the order of 
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100 wppm are Newtonian for shear rates less than about 10^ sec ^ 

(23) . McGary (29) found that POLYOX degrades r^n standing in 
solution as well as under inflijence of turbulence and ultraviolet 
li^t. Degradation in solution also occurs in the presence of 
Vcirious oxidizing agents such as potassium permanganate and hydrogen 
peroxide. The mechanism of this latter degradation is thought to 
be autoxidation initiated by ferrous, cijpric, crprous, and silver 
ions in addition to the agents mentioned above. The quality of tap 
water affects the stability of PEG solutions. Water containing 
chlorine and metallic salts in general yields solutions having poorer 
stability, than distilled water solutions. McGary also confirmed 
that certain alcohols stabilized PEO and reduced degradation. 

Pruit, Rosen, and Crawford (38) observed degradation of WSR-701 
Resin when stored for long periods of time. 
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2. Determination of Drag Coefficients of Spheres 



2.1 Motivation for Experiirents with Spheres . Although much has 
been done in studying drag reduction in pipe flow, conparatively 
little has been reported regarding drag reduction on bluff bodies. 
While the phenomenon is not very well defined for any geoitetry, 
the need for much more e^periirentation with bluff bodies in 
separated boundary layer flow of drag reducing fluids is most 
apparent. 

Pipe flow experiments and rotating disc experirrents usually 
can be carried out without great difficulty; however, experiments 
with bluff bodies present special problems that must be solved 
if meaningful results having general interpretations are to accrue. 
As was mentioned earlier, drag of totally submerged cylinders 
can be greatly affected by L/D. Other complicating problems 
associated with cylinders encountered by McClanahan and Ridgely 
(27) were ventilation behind the strut supporting the test cylinder 
as well as the contribution to drag by the strut itself. 

Ihie fluid containers in which bluff boc^ tests are conducted 
frequently affect the results in some way. For creeping flow, 
very low Pei'Tiolds numbers, the presence of bounding surfaces within 
several diameters from the test body can cause a significantly 
higher value of drag coefficient to result (19,28) . For tests with 
separated laminar boundary layers on bluff bodies, the walls of a 
tunnel cause turbulence in the fluid which may alter the character 
of the laminar boundary layer. This effect has been particularly 
noticeable (15,22,43) for Reynolds numbers in the vicinity of the 
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drag "crisis". According to Torobin and Gauvin (53) it was not 
until the conpletion of the wind tunnel at Goetingen that Vfeiselsberger , 
in 1922, was able bo conduct tests on spheres without being affected 
by turbulence originating at the tunnel walls. 

Atterrpts to surnount this difficulty were mde by observing 
spheres during free fall. While in tunnel tests, the body is sup- 
ported rigidly by a strut and drag force measured directly, in the 
free fa 11 tests terminal speed is determined and the drag force calculated. 
At terminal speed the accelerating force of gravity on a sphere is 
balanced by the retarding force exerted by the fluid plus the buoyant 
force of the fluid: 

(4/3) rr(d/2)^ (d - P.,)g = drag 
s w 

where is the density of the sphere, the fluid density, U the 
terminal velocity and g the acceleration due to gravity. The drag 
coefficient is then: 

^s 2 2 

= (4/3) g d{— - 1)/U = A/Re 

° -"w 

where A is a constant dependent upon sphere parameters and fluid 
density only. It is assuxted that no unbalanced horizontal force acts 
on the sphere. If the sphere is accelerated, indicated by failure to 
reach terminal velocity or deviation from a vertical fall, then the 
alxive expression is invalid. Odar and Hamilton (34) presented an 
analysis for drag on an accelerating sphere in harmonic motion for 
Reynolds numbers below about 60 and later Odar (33) confirmed the 
applicability of the development to freely falling spheres. 
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According to Goldstein (15 ) , Lunnon found that for free 

falling spheres at Reynolds nuinbers greater than about 3 x 10^ 

and less than about 10^ drag coefficients were greater than 

those obtained by wind tunnel tests. Bacon and Reid (2) observed 

5 6 

for Reynolds nimbers between about 3 x 10 and 4 x 10 only a 

gentle decrease in drag coefficient to a minimum value of 0.16 at 

6 

a Reynolds number of 10 instead of the abript decrease to a 

5 

minimum of about 0.08 at a Reynolds nunber 4 x 10 observed by 

others in wind tunnel tests. Lang and Patrick (26) reported 

results of free-fall tests in water consistent with Lunnon 's for 

5 5 

Reynolds numbers up to 2 x 10 . Between 2 and 3 x 10 they reported 
an abrupt decrease in drag coefficient like that obtained by wind 
tunnel testing. For a Reynplds number of 2.5 x 10^, the upper limit 
of their tests, they reported a drag coefficient of 0.15 with an 
apparent decreasing trend. In their tests some of the spheres 
used by Lang and Patrick did not attain terminal velocity. They, 
as did Hayes (20) and Chenard (5) later, corrected their data for 
acceleration by applying a virtual mass coefficient. 

In the present study of drag reducing additives, spheres 
were chosen as the- test shapes and drag determination by free fall 
in a stagnant fluid was chosen as the technique. The choice was 
based on need for sphere ejperimentation especially in the vicinity 
of drag crisis as well as the convenience of having available a 3- 
foot diameter steel cylindrical tank having a height of 6 feet. 

Poly (ethylene oxide)WSR 301 in water was selected as the test fluid 
since, for the range of Reynolds numbers of interest, more results 
are available for comparison. It is apparent that a sphere should 
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attain terminal velocity and a means for determining the value of 
that velocity as well as confirming that acceleration was indeed 
zero should be provided. 

2.2. Methods for Determining Sphere Drag Coefficients . Previous 
investigators using free falling spheres have used various methods 
for velocity determination. Torobin and Gauvin (53) attribute to 
Allen, in 1900, a stereoscopic photogr^hy technique for measuring 
velocities of spheres in fluids. Shakespear, in 1914, made use of 
an electrical tuning system for measuring velocities of small weighted 
celluloid spheres in an air tower, and Bacon and Reid (2) in 1924 
determined velocities of spheres falling from aircraft by using a 
pair of recording theodolites. More recently, incidental to drag- 
reduction studies, Lang and Patrick (26) and Hayes (20) used multiple- 
image strcbe-flash photography for velocity determinations of small 
spheres falling at Reynolds numbers below drag crisis in transparent 
tanks. A. White (64) used high-speed cine photography for larger 
spheres in a transparent tank. In an investigation reported by 
Ruszczycky (41) , Kinnier used a series of coils surrounding a 
cylinder containing the test fluid to sense the change in magnetic 
field as the test shape passed each set of coils. Chenard (5) attached 
to a sphere a very small tether which turned a potentiometer producing 
an output voltage proportional to vertical displacement of the sphere. 
At terminal velocity, the vcltage recorded on an oscillograph yielded 
a straight line plot and the slope of the line corresponded to the 
terminal velocity of the sphere. Thurston and Jones (50,18) tested a 
torpedo shaped body in a 20-foot high tank having nine stations located 
2 feet apart, vertically. Each station had a sensor system allowing 
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detection of passage of the bo<fy. The sensor system was conposed 
of a horizontal transparent slit with a light source shining 
through the slit into the tank and a pair of light sensitive 
resistors mounted level with the slit. These received light 
reflected back from a polished band on the torpedo as it passed 
the station. This system allowed a distance - time relation to be 
obtained and subsequent analysis to obtain velocity. 

A technique employed by Goldstein and Kreid (14) was not 
directly applicable to this study but nevertheless stimulated 
interest in the possibility of altering the technique to fit the 
present task. They msasiired fluid velocity in a duct by measuring 
the doppler shift in laser light scattered from particles moving 
with the fluid. For this study, use of lasers for velocity determin- 
ation was considered by doppler shift rreasurement, as well as by 
eirployment of a laser radar. Both ideas were abandoned because, in 
addition to technological conplexity, the expected speeds of the 
spheres (1 to 9 meters per second) would not likely be determinable 
within the desired accuracy. Other techniques considered but quickly 
abandoned, were acoustic echo-ranging and doppler-shift measurement 
of a reflected acoustic pulse. The use of a photographic technique 
was considered and abandoned since the drop tank was opaque and some 
phenomenal trickery would have been necessary to photograph a reasonable 
field of view from inside the tank. In addition, it was desired that 
a data reduction scheme less cumbersoire than photograph measurements 
could be achieved. It was decided that a technique similar to that 
enployed by Thurston and Jones could be developed and would likely be 
the mast practical of the velocity determination schemes considered. 
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The spheres listed in Table 2.1 were selected in order to yield data 
at Reynolds numbers greater than and less than that at drag crisis as 
well as to provide data within the range of boundary layer transition as 
shown in Fig. 2.1. 

An additional problem which had not been solved was the question 
of how to cause a sphere to achieve teinminal speed in the available 
drop tank. Estimates for the fall distances (S^ required for spheres 
used in this study to attain 0.99 terminal speed shown in Table 2.1 were 
calculated using the ej$)ression developed by Lang and Patrick (26) . 

^0.99 ^1^ 

In this ej^ression is a theoretical virtual mass coefficient given 
as 0.5 for bodies not generating significant wakes and 1.8 when a wake 
was present. The estimates of Table 2.1 assumed = 1.8. 

2.3 Optical Velocimeter . An optical velocin^ter was built as an 
integral unit as shown in Figs. 2.1 and 2.2. Uniformly spaced planes 
of light radiated from two vertical larrp tubes flanking a single 
tube containing a light "receiver". A sphere falling through the 
planes of light reflect lii^t to the receiver throui^ collimating 
slits in an outer shell which surrounded the laitp tubes and the light 
receiver. The light receiver was an acryllic plastic rod contained 
concentrically inside an acryllic plastic tube. Light flashes received 
in the rod were transmitted along the rod to a multiplier phototube 
mounted at the top of the velocimeter. The series of output pulses 
from the phototube as the sphere passed through successive light planes 
could be recorded on magnetic tape or on an oscillograph producing a 
distance - tune relationship for the falling sphere. By dividing the 
distance between two lii^t planes by the time measured between pulses 
a rough mean velocity for that interval was obtained. If terminal 



velocity was attained, this mean velocity is the value for the 
vertical coirponent of the spheres velocity. 

Prior to building the velocimeter described above, a small 
developmental model having only three slits was constructed and 
tested in order to assure feasibility of the velocity determination 
technique, and to develop optimum design parameters. The li^t 
receiver rod was notched at each slit as shown in Fig. 2.4 so as to 
provide for optimum light reception and to allow the acryllic rod 
to act as a light pipe to transmit the maximm amount of light to 
the phototube. It was e^$>erimentally determined that variation in 
the angle of the face of the notch produced variation in the photo- 
tube output. A "V" cut having an included angle of about 40° - 60° 
and a depth of cut approximately 1/6 of the rod diameter gave slightly 
greater output than other notches. Polishing the notch slightly 
diminished the output. Turning the rod so as to position the slit 
to the side opposite the incoming reflected light produced approxi- 
mately 30% increase in voltage output of the phototube. It was also 
determined that without the notch, no significant light reflection 
was received at the phototube. Misalignment of the notches behind 
the collimating slits by as little as 0.5 mm caused drastic reduction 
in the phototi±>e output-pulse amplitude. This caused serious concern 
because the signal-to-noise ratio was low (about 15 db) even under 
optimum velocimeter operating conditions. 

Light pipes are optical waveguides in that because of multiple 
internal reflection light trapped inside is transmitted along the 
"pipe" even though slight bends may occur. Two transparent media 
having a smooth interface and having different indices of refraction 
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will allow internal reflection. If a light ray incident on the 
interface in the itiedium having the hi^er refractive index, n ' , 
has an angle of incidence greater than the critical angle, (J>^, 
given by: 

^ • -1 , . 

= sin n/n 

where n is the refractive index of the second rredium, then the 
ray will, in principle, be totally reflected. In practice, 
iiiperf actions in the interface due to defects or oontamination 
by dust and grease cause some light to be transmitted or leaked 
through the interface. For a short light pipe or even a long pipe 
transmitting very intense light this might not be a problem. Since 
the light reflections in the velocimeter would be of low intensity 
It seemed necessary to optimize the light receiver surface. The 
possibility of silvering the rod that would, in the final design, 
likely be approximately 6 feet in length, was considered. It was 
hoped that this would prevent contamination of the plastic surface 
and "help" the reflectivity somewhat. A commercial firm was consulted 
and because of the nature of the silvering process the acryllic plastic 
rod could be silvered only by very special treatment V\Aiich, it was 
estimated, would cost on the order of $9000. Because of the cost 
involved and because it had been reported by Siegmund (47) that silvering 
could actually result in greater li^t leakage, the idea of a metallic 
coating was abandoned. Siegmund also reported that the best coatings 
for light pipes were dielectric substances having refractive indices 
less than that of the pipe but because of the e^^cted small advantage 
no coating was specified for the velocimeter light receiver. For the 
velocimeter a low value of <t>^ was desired for maximum light reception 
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or "trapping" and transinission. Since = sin~^ n/n' the ratio 
of n/n' had to be minimized. Since the refractive index of the 
plastic was 1.5, that for water was 1.33, and that for air was 
1.0, it was decided to encase the light receiver inside a water- 
tight acryllic-plastic tube viiich would allow the light pipe to 
function with a plastic-air interface and would give nore favorable 
internal reflection than would a plastic-water interface. The rod 
would have considerably more than a few wavelengths separating it 
from the tube and thus no significant light leakage was ejected. 

The light source for the developmental velocimeter model was 
a pair of fluorescent tubes, one in each lairp cylinder. Immediately 
\;pon atteirpting to operate the device it was apparent that the electro- 
magnetic noise generated and radiated by the fluorescent process was 
excessive and conpletely obliterated any reflection pulse that might 
have been present in the phototube output. It was determined that 
a d.c. lamp would allow satisfactory operation if an adequately high 
intensity lamp could be found. The highest intensity d.c. lamps 
available with globe size of one-inch diameter or less was a 12-volt 
automobile instrument lamp, number 1073, \diich had a rating of 32 
candlepower. The one-inch globe diameter was dictated because the 
lamp cylinders had been designed for a one-inch diameter fluorescent 
tube and, in addition, small lamp-cylinder size was desirable for 
good collimation with reasonably small overall diameter. Associated 
with the choice of light source was the problem of noise in the photo- 
tube output even with a d.c. light source. In order to obtain 
satisfactory velocimeter operation it was determined that the use 
of a low-pass filter connected to the photcit^±e's cathode follower 
output was necessary. No satisfactory explanation of this is available 
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but apparently the very low light intensities reflected into the 
light receiver were inadequate to overcoire dark current noise and 
other electronic noise in multiplier-phototube operation and the 
velociiteter pulses were distinguishable only by filtering. It was 
observed that if a sufficiently intense light source were beamed 
into the slits, "normal" phototube operation (without filtering) 
was possible, confirming, to some degree, the conjecture given 
above. Tests were condqcted in which laitp intensity was varied by 
varying sipply voltage, high voltage supply to the phototube was 
varied, and filter cutoff frequency was varied. For the model 
velocimeter, it was determined that the conbination of parameters 
giving best results was: rated voltage on each of the lamps (12 

volts) , 1550 volts high voltage supply and, low pass 35 Hz and 
below in the filter. In other tests it was noted that the phototube 
output-pulse arrplitude decreased very rapidly with separation between 
velocimeter and sphere. This decrease was found to be approximately 
proportional to the square of the sum of the distance between the 
light source and the sphere and between the sphere and the light 
receiver. 

A Dumont 6291 multiplier phototube was used in the model 
development tests but it was determined that, although a Dumont 6292 
had nearly the same characteristics, sli^tly irtproved velocimeter 
performance resulted by using the 6292 phototube. This latter tube 
is a 10 - stage multiplier phototube with a flat end - window type 
photocathode having its maximum response at about 4400 Angstrcms. 

The response of this multiplier phototube was such that incandescent 
laitps with tungsten filaments provided a very good match with the chosen 
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multiplier phototube. De Vos (6) reported that the emissivity of 
a tungsten ribbon decreased very slightly in the wavelength range 
between 0.35 u and 0.55 m / which is the range of greatest sensitivity 
of the Dumont 6292 phototube. In this range, according to Pluck's 
radiation formula, radiated energy increases with increasing wave- 
length and thus the laxrps' spectral emission matched the spectral 
response of the phototube reasonably well. The Duiront 6292 multiplier 
phototube has an adjustable focusing shield allowing for optimum 
collection of photoelectrons. It has a riimimimi window diameter of 
1.5 inches and a window index of refraction of 1.5 v\iiich is the sane 
as that of the acryllic plastic rod from vdiich light is transmitted 
to the phototube. 

During testing of the model velocimeter the duration of the pulses 
produced as a sphere passed a light plane was approximately 10 to 12 
msec. This is a much longer pulse duration than e^^cted from slit 
geometry considerations and was probably caused by the use of the 
narrow-band filter at the photomultiplier output which led to longer 
circuit rise and decay times than would have resulted had no filter 
been used. Since it was desirable that each pulse be distinguishable, 
the required separation of light planes for a sphere falling at 8.6 n\/sec 
was approximately 0.1 iteter. Twelve light planes uniformly spaced at 
0.108 m were specified for the final velocimeter as a coirpromise of the 
following : 

Ij Provide accuracy for measurement of time between pulses by 
large separation of the light planes. 

2) Provide a maximum number of pulses to allow more samples for 
distance-time recording, and thus better resolution of velocity 
changes. 
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3) The velocinneter should not extend to the bottom of 
the tank since a sphere rebounding from tank bottom 
could cause extraneous signals which could interfere 
with analysis of valid pulses. 

4) Light receiver rod should be of minimum length to avoid 
excessive light leakage and attenuation. 

The optical velocimeter as shewn in Figs. 2.1 and 2.2 had 
twelve 12 volt d.c. lartps in each lairp tube in 3 grovps wired in 
series, each group consisting of four latrps in parallel. This 
allowed each lairp tube to be operated with its own power siprply at 
36 volts and 6.8 arrperes. The laitps in each tube were mounted on 
a plastic strip as shown in Fig. 2.5 so that each lanp socket could 
be moved in its mounting block. This provided capability of positioning 
each lamp filament for uniform spacing or whatever filament spacing 
was required. In addition the plastic strip could be adjusted vertically 
for optimum position of the entire lamp asserrrbly. Provision was made 
for a flow of cooling air within the lamp tubes to allow extended 
operation of the laitps if necessary with minimum danger of damaging the 
acryllic plastic tube and the insulation on the electrical leads to 
the laitps. 

Optimum operation of the velocimeter was obtained by setting lamp 
tube voltage at 36 volts, multiplier phototube high voltage at 1050 
volts, and the SKL Model 302 variable electronic filter set to pass 
frequencies of 120 Hz and below. The velocimeter outer casing was 
coated inside and outside with a light absorbing and diffusing velvet 
black paint. The aluminum lamp cylinders containing the lamp tubes as 
well as the aluminum cylinder containing the light-receiver assembly 
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was given the saire finish. This was to prevent extraneous light 
reflections from interfering with proper velocimeter operation. 

The inultiplier phototube was mounted inside an aluminum cylinder 
which also contained the cathode follower circuit for proper 
matching of the output with other electronic equipment used in 
conjunction with the system. A light-tight stainless-steel flange 
was designed to attach to the light-receiver cylinder and thread 
into the photomultiplier's aluminum cylindrical housing. 

2.4 Sphere Ejector. In order that the test spheres achieve terminal 
velocity within the 6-foot drop tank, it was necessary to provide the 
spheres with an initial velocity. The optical velocimeter was to be 
positioned with its first light plane 6 to 8 inches beneath the fluid 
surface and it was desirable but perhaps not realizable to have the 
sphere at its terminal velocity as it passed the first light plane. 

An air-driven sphere-ejector system shown in Figs. 2.2 and 2.6 was 
built which catapulted a sphere downward, accelerating it from zero 
velocity just beneath the surface to near terminal velocity as the 
driving piston cortpleted its downward travel. The ejector consisted 
of a cylinder and piston, sphere holder, accumulator, pressure regulator, 
and qmck-action solenoid-operated valve. To the shaft of the piston 
was attached a short nylon rod which contacted the sphere. Sphere 
holders or "claws" of three sizes to accommodate the varying size 
spheres were constructed. Each had an aluminum base and three shaped 
phosphor bronze spring clips which held the sphere in position until 
the ejector shaft pushed the sphere from the holder. The sphere was 
ejected from a point 5 to 10 cm beneath the fluid surface so that 
"bubble capture" would be no problem. 



45 



2.5 Data Recording and Processing . Calculation of accurate drag 



coefficients and Reynolds nunibers required maximum accuracy in 
measuring the time between recorded pulses and required that the 
distance between light planes be known accurately. The tine 
between pxiLses was determined by recording the velocirteter output 
during the fall of the sphere at 37.5 inches-per-second (ips) on 
a Ptecision Instruments Model PI-6200 Instrumentation Tape Recorder 
operated in FM Mode. Ihe recording was then played back at 3.75 
ips and the tape recorder output was connected to a Brush Instrument 
Model BL-202 Oscillogr^h through a d.c. anplifier. The oscillograph 
was operated at 125 mm/seoond allowing the pulses to be re-recorded 
with a greatly expanded scale so that an accurate measurement of 
distance between pulses on the oscillogram could be made, which 
corresponded to a measurement of time. With the playback at 1/10 
real-time speed and the oscillograph at 125 mm per sec, a scale of 
0.0008 sec/mm was established. However, it was recognized that 
mechanical recorders mii^t not always run exactly "on speed". In 
order to provide a check for any speed variation in the recording- 
playback link, a sinusoidal time-referenoe signal from a General 
Radio 1309A Oscillator was recorded on the second channel of the 
magnetic-tape recorder simultaneously with recording of the velocimeter 
output on the first channel. This reference signal was then recorded 
on the second channel of the oscillograph viiile the velocimeter pulses 
were recorded on the first channel. The resulting oscillogram was 
coitpletely independent of any recording or playback speed variations. 

It was necessary, however, to insure that the known reference signal 
was set properly for recording and that the period of the reference 
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signal on the oscillogram be measured for each interval between 
pulses. Since the recording equipment was high quality it was 
decided to sacrifice some flexibility in the interest of carputational 
convenience and assume that the recording-playback speed during the 
fall of a sphere could be approximated with sufficient accuracy by 
an average speed for that oscillogram. Thus, instead of correcting 
the time for each interval transit with its own unique correction 
factor, a time correction factor common to the entire oscillogrcim 
was applied. As a further conputational convenience, a unique time 
reference signal was assigned for each sphere used in the investigation. 
This allowed identification of a particular sphere for an oscillogram 
by measuring the period of the reference signal. 

Since the carputations in this study would soon become unmanageable, 
a program for handling the sphere data by an IBM 360/67 digital conputer 
was written and is listed in Appendix A. Ideally the data would have 
been processed by playing the magnetic-tape recording into the analog- 
cortputer section of a hybrid-corrputer system such as the SDS 9300/COiyDDR 
5000 system. The analog conputer 's logic could have been enployed to 
count time between pulses and pass the information into the digital 
conputer for real time conputation of velocities, drag coefficients, 
and Reynolds numbers. The digital conputer program given in Appendix A 
was originally written with hybrid operation in mind but the analog logic 
was not conpleted, in time. 

The digital computer program inputs for each oscillogram correspond- 
ing to one fall of a sphere were; 1) eleven pulse interval measure- 
ments in millimeters, 2) eleven reference signal periods in millimeters, 
3) fluid density in kg/m^, 4) shear viscosity of the fluid in poise 
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corrected for tenperature, 5) first run number for the batch 
of oscillograms being processed, and 6) any runs suspected of 
containing errors and which should not be used in conputing 
averages of drag coefficients and Peynolds numbers. 

The program contained information for each sphere and upon 
processing the reference signal information the information for 
the proper sphere was selected and used with a calculated average 
velocity in oortputing the drag coefficients and Reynolds number for 
that sphere for that run. The program also contained the interval 
distances between each pair of light planes vAiich had been calculated 
from numerous measurements to assure precision in this critical part 
of the velocity determination. VJhen several runs for a sphere were 
processed together, averages for drag coefficients and Reynolds 
numbers for that sphere were also calculated. As an aid to data 
analysis the punched data cards were returned to the conputer with 
a program which caused a velocity history, plot of velocity for each 
velocnreter interval, to be printed. This allowed visual assessment 
of whether or not a sphere attained terminal velocity during its fall 
as well as detection of any unusual velocity behavior during the fall. 
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3. Ej^rimental Procedure 



3.1 Apparatus . The optical velocirteter was installed in 

the drop tank and its associated equipment connected as shown in 
Figs. 2.6 and 3.1. The sphere ejector was mounted on a heavy 
steel cantilever support above the tank as shown in Fig. 2.2. 

The solenoid which activated the air valve of the sphere ejector 
was actuated by a multiple contact switch which also connected 
the time reference signal to the magnetic tape recorder. This 
served as a convenience for finding the pulses on the magnetic 
tape. A portable air coirpressor capable of delivering small air- 
flow rates at pressures up to approximately 145 psi was used as 
a source for the sphere ejector. The regulator and accumulator 
allowed setting any desired pressure for sphere ejection. Since 
the phosphor bronze "claws" held the spheres quite firmly a minimum 
"break out" pressure was necessary to release the sphere from the 
holder. Thus it was not possible to utilize pressures much less 
than about 10-15 psi which, in some cases, led to initial sphere 
velocities considerably in excess of tdie desired terminal velocity. 

The distance between light planes was meas^ared by suspending 
a tethered 3" aluminum sphere over a pulley into the drop tank 
filled with water. The pulley was placed so that the sphere could 
be lowered into the tank along essentially the same path it would 
take were it catapulted from the ejector. The tether was marked and 
mowd along a rteter stick so that as the sphere was lowered into 
the tank the depth of the sphere could be neasured with reasonable 
accuracy. As the sphere passed through a light plane the velocimeter 
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output pulse was displayed on an oscilloscope allowing the sphere 
to be positioned for each light plane at that point at which iteximum 
light was reflected into the light receiver. The depth of the sphere 
corresponding to the "center" of each light plane was determined 
and the difference in depth fron one plane to the next was taken as 
the interval distance. Several sets of determinations were made in 
order to obtain reasonable assurance that these critical parameters 
given in Table 3.1 were properly precise. Sphere ejector performance 
was evaluated in order to ascertain proper ejector pressure for the 
spheres to attain terminal velocity. Sphere number one was ejected 
in water at various pressures and the mean velocities within the first 
interval of the velocimeter were calculated as shown in Fig. 3.2. 

3.2 Procedures 

3.2.1. General . For the drag-coefficient determinations the 
ejector pressure was set to the appropriate pressure which would 
cause the sphere to quickly attain its "water" terminal velocity. 
Pressure adjustments could be made fron one set of rms to the next 
if necessary. 

3.2.2. Tests in Water . To develop confidence in the techniques 
enployed in this stuc^, drag coefficients of the spheres falling in 
water were determined and comparison with accepted values was made. 
The weight in air and diameter of each spheie was determined, the 
density was calculated, and the results, shewn in Table 2.1 were 
incorporated in the digital computer program. Before beginning a 
set of runs all electronics equipment was allowed ample warm-up time, 
but high voltage was not applied to the multiplier phototi±)e. A 
sphere was loaded into the ejector and the proper pressure was set 
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as described above. The high voltage to the multiplier phototube 
and the laups were turned on. A few seconds delay was allowed, 
to insure the decay of transients in the cathode follower circuit 
due to the application of high voltage and light. The tape recorder, 
operated in FM mode, was placed in "record" and allowed to 
accelerate to "on speed" at 37.5 ips. The elector switch was 
actuated causing the sphere to be ejected and the reference signal 
to be connected to the tape recorder. The recorder was stopped 
and lairps and high voltage turned off. Upon corpletion of several 
runs the fluid tenperature was measured and recorded. The tape was 
then rewound and played back at 3.75 ips, the output recorded on an 
oscillograph, and the data analysis conpleted as described previously. 

3 2 3 Tests a n Poly (ethylene oxide) Solutions .sphere drag 
coefficient determination in very dilute PBO solutions was carried 
out in the sane manner as that for water except, of course, fluid 

preparation provided a greater challenge . 

3.2.4. Mixing S olutions . Although POLYOX resin was water 
soluble, it coiald not be dissolved in water in the same manner as 
an inorganic salt. Because of the tendency for the surface of a 
particle to becote soft and sticky in the presence of water, if the 
resin particles were not dispersed prior to contact with water they 
tended to form a gel-like agglomeration which could only be dissolved 
after a very extended period of time. In order to dissolve PEO resins 
in water. Union Carbide Chemicals Company (59) recoranended pre- 
dispersion be acconplished either by mechanical means or by taking 
advantage of the insolubility of PEO resins in various liquids. In 
the latter case boiling water could even be used as the dispersion 



51 



carrier due to the inverse tenperature - solubility property of PEO. 

Of more practical use, however, was the use of a polypropylene glycol 
for the dispersal agent as was done by Forester and Francis (10) 
during friction reduction tests for the David Taylor IVbdel Basin 
and by McClanahan and Ridgely (27) . The polypropylene glycol used 
in their tests was of molecular weight about 400 and was marketed 
by Dow Chemical Company as POLYGLCOL P400. PEO was corrpletely 
insoluble in POLYGLYCOL but the POLYGLYCOL was miscible with water 
allowing the resulting dispersion to be injected into water so that 
the PEO could be readily dissolved. Of significance is the observation 
made by Forester and Francis that the dispersion carrier produced no 
drag reduction when injected alone into water. Thus, any drag reduction 
observed was due to the PEO and not the polypropylene glycol. For 
this investigation PEO concentration of 1000 , 200 , and 100 wppm were 
to be tested. A slurry of 1.081 Kg of PEO suspended into 2 Kg of 
POLYGLYCOL was prepared. This slurry, when dispersed into the drcp 
tank containing 1081 Kg of fluid when filled with tap water to with- 
in 0.1 ra from the top, yielded an aqueoios solution of 1000 wppm 
poly (ethylene oxide). After the dispersion was corrpletely errptied 
into the tank and the tank was filled to the proper level the solution 
was stirred gently for 18 hours with two Cole-Parmer mixers operated 
at very slew speed to insure corrplete solution of the resin without 
undue polymer degrading turbulence and to provide some assurance of 
homogeneity of the solution. For tests with 200 v^pm a slurry of 0.216 
Kg PEO in POLYGLYCOL was prepared and etrptied slowly into the tank while 
it was being filled with tap water. Gentle stirring by both mixers was 
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carried out during the entire time the slurry was enptying into the 
tank as well as for about two hours afterward when the PBO resin 
appeared to be conpletely dissolved. The 100 wppm PBO was prepared 
by draining the 200 wppm solution from the tank to a fluid level 0.936 m 
below the tank top, refilling with tap water, and stirring. 

It has been reported by McGary as well as other investigators 
that several ions often present in tap water cause degradation in 
POLYOX. A typical analysis of the water used is given in Table 3.2 
in order that some assurance be given that no significant degradation 
due to chemical action affecting drag reduction would result. 

After the solution was prepeured sphere drag coefficients were 
determined as was done for water. For each concentration sarnples 
of the fluid were taken at several levels in the tank and capillary 
tube viscosity determinations were made. This allowed obtaining 
infoimation as to the actual homogeneity of the fluid as well as to 
determining intrinsic viscosity of the POLYOX resin. The molecular 
weight of the PBO resin used was estimated using the ej<pression given 
by Shin (46) relating intrinsic viscosity to molecular weight for 
25°C: 

r , T ir,-4 „,0.78 

[r] = 1.03 X 10 M 
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4. Results and Discussion 



4.1. Drag Coefficients in Water and Dilute Poly (ethylene oxide) 
Solutions . The results of the sphere drag-coefficient determinations 
in water, 1000 wppm PEO, 200 vppm PEG, and 100 wppm PEO are presented 
in Fig. 4.1, 4.2, 4.3 and 4.4 respectively. Since no significant data 
display iitprovement would have resulted by pretending that the 
solution shear viscosity obtained by capillary-tube visconetry had 
meaning, drag coefficients were plotted with water Reynolds numbers 
instead of solution Reynolds numbers. This was done as a matter of 
convenience for coitparing these data with the results obtained by 
other investigations. Except for the 1000 wppm solution, the resulting 
solution Reynolds numbers would have differed from the water Reynolds 
numbers by a barely perceptible ainount. The use of the generalized 
Reynolds numbers for power law fluids suggested by Dodge and Metzner 
(8) would not have resulted in an urproved data display. 

The wind-tunnel data for spheres attributed to Vteiselsburger 
as presented by Schlicting (43) are shewn as a solid curve on all 
the graphs. 

It is significant that the sudden drop in drag coefficient 
associated with the drag crisis observed for spheres moving in air 
and water was also observed at about the same Reynolds number for all 
PEO solutions investigated. However, the values of drag coefficients 
for Reynolds number less than the critical value, as well as those 
slightly greater, were considerably lower than for water and air. 

Data at the lower limit of Reynolds numbers for this investigation 
coincided with results reported by Lang and Patrick at the upper 
Reynolds nurtber limit of their investigation. 
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Another significant observation is that vAien the Reynolds 
number is increased above the critical value/ the drag coefficient 
returns precipitously to nearly the same value as for water. This 
feature could not be further defined since it occurred at the upper 
limit of Reynolds nuirbers for this investigation, about 8 x 10^. 

In the region above the critical Reynolds number a very slight PBO 
concentration dependence was observed with the 1000 wppm producing 
greatest drag reduction and 100 wppm producing least drag reduction. 

Sphere size did not affect drag reduction if the water Reynolds 
nunber was equal. It is seen that the 2 1/2 inch brass sphere and 
4 inch aluminum sphere had approximately the same drag coefficient 
at a Reynolds number of about 6 x 10^ in 1000 wppm PEO. 

4.2. Sphere Behavior During Fall . During the fall of seme spheres , 
deviation from the vertical path, either curvature or occassional 
spiralling was observed. For a given sphere in a particular fluid 
the deviation seerted to be the same from one fall to the next. The 
maximum eccentricity of the spheres was 0.2% and this was not a 
contributing factor. The sphere holder on the ejector was readjusted 
and the same deviation occurred. All spheres fell more smoothly and 
had more nearly vertical paths in 1000 vppm PEO. The velocity cor- 
rection would have been only about 1.03 for a sphere moving along a 
straight path from the tank center at the top to the tank side at 
the bottom. However, those runs in vhich the sphere deviated widely 
from the vertical path were excluded from the data presented here. 

The brass spheres, almost without exception, fell in vertical paths. 
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As had been anticipated (see Fig. 2.1) the 3 inch aluminum 
sphere, (number 8), yielded two values of drag coefficients. 

Other spheres had been ejq^ected to yield multiple results but 
sphere 8 was the only one tested specifically for this effect. 

The ipper value of drag coefficient (see Figs. 4.1 and 4.2) was 
obtained by setting the sphere ejector for minimum ejection speed 
and the lower values by selecting a higher ejection speed. A 
phenomenon of interest was observed in the speed histories. In 
the first cases the sphere decelerated to its terminal speed. 

For the lower drag coefficients, the sphere accelerated and 
settled to its terminal speed with no observable oscillation. 
Apparently, in the second case, the ejection speed had been 
sufficiently high to sJcip over the lower terminal speed. It is 
likely that a sphere falling from rest could not have made this 
"junp", althou^ no supporting evidence for this was obtained. 

This phenonenon is not related to the dual terminal velocity 
behavior for sate spheres reported by Sanders (42) since the 
present phenomenon was predictable on the basis of the sh^)e of 
the drag coefficient curve tdiereas that reported by Sanders 
apparently resulted from another mechanism characteristic of the 
fluid. 

The paucity of data presented for 200 wppm and 100 wppm PEO 
resulted from the exclusion of so many drops in which the spheres 
did not decelerate to terminal speed after having been ejected at 
a much too great ejection speed. This was not apparent until after 
data analysis begcin and the plots of time history for each drop were 
obtained from the computer. Had these histories been available sooner 
it might have been feasible to repeat many drops with mare appropriate 
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ejection speeds. Nevertheless, even with so few data, the general 
trend indicated by the 1000 wppm data, can be observed in the 200 
wppm and 100 wppm data. 

4.3. Degradation of Poly (ethylene oxide) . No degradation resulting 
from the spheres falling in the PEO solutions could be determined. 
Intrinsic viscosity determinations by capillary tube yLscotetry, 
before and after a days run produced no measurable change in the 
value of intrinsic viscosity from the initially nEasured value of 
8.6 j;_ 1*0 deciliters/gram at 25°C and thus probably no change in 
molecular weight. In addition, the shear viscosity of 200 v\p>pm 
solution samples changed from 1.13 centipoise at 25°C within a 
few hours after mixing to 1.12 centipoise at 25°C ten days later. 

This indicated that the likelihood of degradation due to ageing 
during the sphere drops conducted within 24 hours after mixing was 
very small. 

The expression given by Shin (46) for PEO molecular weight 
and intrinsic viscosity at 25°C was used to calcxilate the molecular 
weight of the PEO WSR 301. The resulting value of 2.0 x 10 is 
approximately half that given by the manufacturer. The resin used 
had been stored in its unsealed paper carton for approximately 2 1/2 
years prior to this investigation. Polymer resin degradation cannot 
be proven on the basis of this determination alone since no molecular- 
weight determination for the new resin was available. Since, Pruit, 
Rosen, and Crawford (38) obtained a value of 1.57 x 10^ for WSR 301 
there may be no problem here at all. 
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5 . Conclusions 



5.1 Drag Coefficients of Spheres . The effect of the addition of 
poly (ethylene oxide) to water on the drag of free-falling spheres 
is best discussed with the assistance of a plot of drag coefficients 
vs water Reynolds number (Fig. 5.1). The lines represent a "best" 
estimate of the shape of the drag coefficient curve as visually 
fitted to all available data. These estimates were based on the 
results reported by Lang and Patrick (26) , Hayes (20) and Chenard (5) 
for Reynolds number less tlian 1.5 x 10^ and on the results of this 
investigation at greater Reynolds numbers. 

From results reported by previous investigators it is seen 

4 

that, as Reynolds number increases above about 10 , drag enhancement 
changes to drag reduction which increases with increasing Reynolds 
number. Furthermore, drag reduction is increased by decreasing 
concentration from 1000 wppm to 100 V\^pm. 

As the critical Reynolds number is approached from below, the 

5 

concentration dependence seems to become less until about 2 x 10 
it is not measurable. Admittedly, this latter conclusion is based on 
a minimum of data but, because of the very great difficulty in obtaining 
free-fall data at the critical Reynolds number, it is not likely that 
the curve will be better defined unless ej$>erimental technique superior 
to those used in the past are developed. Nevertheless, the results 
of this investigation show that at just above 2 x 10^ there is no 
measurable concentration dependence for sphere drag coefficients in 
1000 wppm and 200 wppm PBO solutions. The estimated curve from this 
point until Reynolds numbers above the critical value are reached is 
a guess since there is no data, except for the end points, to support 
the curve shape estimate during the drag crisis. 
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In the Reynolds nurrfoer region between 3 x 10^ and about 7 x 10^, 
there is a very sirell polyrter concentration dependence with 1000 wppn 
PBO solutions giving greatest drag reduction and 100 wppm solutions 
the least. 

An abrupt rise occurs above 7 x 10’ which suggests a behavior 

similar to that of cylinders as measured by Roshko (40) and shown 

in Fig. 5.2. This lends further support to the speculation by 

Sanders (42) as to the applicability of RosWto's laminar separation 

bubble to flows about spheres. If this speculation has any basis, 

then the drag reduction ei^Jerienced in the Reynolds number region 
5 5 

between 3 x 10 and 8 x 10 , called the supercritical region by 
Roshko, IS characterized by the effect of the polymer on the laminar 
separation bubble. This explanation of pol^Ttier effect in the super- 
critical region is consistent with the mechanism proposed by Sanders 
(42) for subcritical Reynolds numbers. In both cases the polymer 
effect IS restricted to altering the character of the laminar boundary 
layer . 

5 , , . 

At a Reynolds nunt>er of about 8 x 10 , drag coefficients in all 
fluids tested in this investigation (including water) were of nearly 
the same value. This apparent equality of drag coefficients may be 
indicative of turbulent boundary layer separation and, further, 
indicative that the turbulent boundary layer is unaffected by polymer 
concentration . 

As indicated in Fig. 5.1 the use of the virtual mass coefficient, 
K,= 1.8, for separated laminar boundary layer flows past spheres, as 
was done by Lang and Patrick (26' , Hayes (20) , and Chenard (5) , was 
corroborated. No such support was given, however, for the use of 
K,= 0.5 for separated turbulent boundary layer flows as used by (Thenard. 
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Sphere size had no effect on drag reduction in PBO solutions 

5 

at a water Reynolds number of 6 x 10 . 

5 . 2 Drag Coefficients of Free-falling Spheres in Water . Reference 
to Fig. 4.1 shows that below drag crisis the drag coefficients 
obtained in water agree more closely with the results of Lunnon's 
investigations of free-falling spheres (see Goldstein 15 ) than they 
do with l^iselsberger ' s wind tunnel data. On the other hand, beyond 
drag crisis the drag coefficients are still slightly greater than 
those reported by Vfeiselsberger but much less than those reported 
by Bacon and Reid ( 2) . It has long been thought that, for free- 
falling spheres, the Lunnon curve could be extrapolated across the 
drag crisis region to join the Bacon and Reid curve. The results 

of the present investigation are inconsistent with this extrapolation. 
It spears that either the ejected spheres of this investigation did 
not behave strictly as free-falling spheres or that the Bacon and 
Reid results cannot be valid for free-falling spheres in the super- 
critical region. It may be that sane slight tiurbulence in the drcp 
tank may have altered the fall in some way but it is obvious that 
the spheres falling in the atmosphere as investigated by Bacon and 
Reid must have been influenced in sa:tie similar way. Thus, it seems 
that beyond drag crisis good wind tunnel data, such as that of 
Weiselsberger , better describes the drag of free-falling spheres 
than does the Bacon and Reid results. 

5.3 Polyner Degradation . The^^^(ethylene oxide) solutions did not 
measurably degrade during the sphere drops nor was degradation due to 
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aging of a 200 wppm solution observed within 10 days after the 
solution was itdxed, 

5.4 Coiments on the Ej^ruTent . The technique of ejecting spheres 
vertically at or near terminal speed and measuring the vertical 
velocity corrponent provided an effective means of extending the 
range of ej^riments in PEO solutions. The limiting factors were 
significant deviation from a vertical path by some spheres during 
soitB drops and the laborious data reduction involved in reading 
oscillograms and manually inserting digital data into the conputer. 
Of mere than 300 drops data analysis showed that only about 150 
could be used with certainty due to erroneous selection of ejection 
speed. Nevertheless, if one is persistent, this technique can 
produce iteaningful results for spheres in fluid flows at and above 
the critical Reynolds number of water. 
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WEIGHTS OF POLYOX RESINS 
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0.05083 + 0.00010 



TABLE 3,1 



DISTANCES BETWEEN LIGHT PLANES OF OPTICAL VEIOCIMETER 



Ihterval, I 



Bounded by light planes Distance, S(I),m 



1 


1 


- 


2 


0.1078 


+ 0.0064 


2 


2 


- 


3 


0.1075 


+ 0.0049 


3 


3 


- 


4 


0.1092 


+ 0.0043 


4 


4 


- 


5 


0.1103 


+ 0.0035 


5 


5 


- 


6 


0.1101 


+ 0.0029 


6 


6 


- 


7 


0.1077 


+ 0.0019 


7 


7 


- 


8 


0.1101 


+ 0.0027 


8 


8 


- 


9 


0.1073 


+ 0.0038 


9 


9 


- 


10 


0.1109 


+ 0.0034 


10 


10 


- 


11 


0.1102 


+ 0.0027 


11 


11 


- 


12 


0.1093 


+ 0.0032 



TABLE 3.2 



TYPICAL GRAVIMETRIC ANALYSIS OF WATER 
SUPPLIED BY IDNTEREY CITY WATER COflPANY 



Substance 


PPM 


Substance 


PPM 


Silica 


21.1 


K 


2.0 


Iron 


0.1 


00^ 


99 


Mn 


0 


Cl^ 


13.0 


A1 


0 


SO. 


28.0 


FI 


0.25 


Total dissolved 




CO_ 


3.1 


solids 


165 


Ca^ 


26.4 


Total CaCo^ 


105 


Mg 


9.5 


NonCarbon 




Na 


11.9 


hardness 


24 



PO . 0 . 05 

7.7 pH 
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Pseudoplastic 




Figure 1.1 Shear Flov Curves 




Figure 1.2 Friction Factor Design Chart for Newtonian and Non- 
Newtonian Fluids. 
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Figure 1.3 Friction Reduction in Pipes for Some Additives. 




65 





Water Reynolds Nunber 



Figure 1.6 Torque Reduction for a Rotating Disk. 
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Figure 1.7 Drag Ccefficients for Circular Cylinders (Wsiselsberger) . 




Figiire 1.8 Drag Coefficients for Spheres (Weiselsberger) . 
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Figure 1.9 Trend of Drag Coefficients for Circular Cylinders in 
Dilute Aqueous PEO Solutions. 
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1.10 Trend of Drag Coefficients for Spheres in 
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Reynolc3s Number 

Figure 2.1 Predicted Re Flanges for Spheres of Various Sizes and Materials. 



QUICK-ACTION SOLENOID VALVE 




figure 2.2 View of Top of Drop Tank with ^here Ejector 
and Optical Velociiteter Installed. 
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SPHERE HOLDER 
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Figure 2.3 Disassei±)led Optical Velocinieter 



Aluminum outer Light receiver 




Figure 2.4 Diagram of path of li^t rays into light reoeiver. 
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Figiire 2.6 View of Apparatus Enployed in Sphere Drag 
r-feasureinents . 
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Ejection Velocity, iti/sec 




Ejector Air Pressure, PSI 



Figure 3.2 



^here Ejector Performance with 4 Inch Brass Sphere. 
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Figure 5.2 Measured Circular Cylinder Drag Coefficients. 
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